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ABSTRACT 



o 

u 

o 



A method is proposed for deriving the dispersion measure (related to the hne- 
of-sight integrated free electron density) to high-redshift, powerful synchrotron radio 
sources by analysis of the temporal/spectral properties of the recorded signals. This 
would allow direct measurements of the distributed density of ionized baryons along 
the line-of-sights to the sources. 
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1 INTRODUCTION 

Highly dispersed radio pulses of likely extragalactic origin 
have been detected at GHz frequencies - one with a disper- 
sion measure DM= J d£ rie — 375 pc(cm)~^ (Lorimer et al. 
2007) and another with DM= 746 pc(cm)"^ (Keane et al. 
2012). Here, the integral is from the earth to the source and 
Ue is the free electron density. Unfortunately the pulses do 
not appear to repeat and the distances to their sources are 
unknown. Here, we propose that the dispersion measures to 
high-redshift incoherent synchrotron radio sources can be 
extracted by analysis of the amplitude cross-correlation of 
the radio waves received at two nearby frequencies. This 
method utilizes the broadband nature of the synchrotron 
radiation. 

Section 2 discusses the theory first without intervening 
plasma and then including it. Also, the influence of inter- 
stellar/intergalactic scintillations are considered. Section 3 
gives our conclusions. 



2 THEORY 

The measured electric field on the ground from a distant 
optically-thin incoherent synchotron source can be approxi- 
mated as 



(1) 



where the different terms represent the contributions from 
the individual electron or positron in the source, £j is the 
complex field amplitude, and S is the delta function. The 
synchrotron radiation from a highly relativistic electron with 



Lorentz factor 7 ^ 1 is directed in a narrow beam of angular 
width ~ 1/7 in the direction of the electron's instantaneous 
velocity (e.g., Pacholczyk 1970). This narrow beam sweeps 
across the observer in a time interval dt ~ {■y^ujc)~^ which 
is much less than the gyration period of the electron, 2Tr/uJc- 
Here, luc = eH / (m^'yc) is the relativistic cyclotron frequency 
with e and the charge and rest mass of the electron and 
H the magnetic field strength in Gauss. Thus the radiation 
at a point on the ground is in narrow pulses each of width 
~ dt which is the basis for using delta functions in equation 
(1). Consequently the frequency spectrum of the radiation 
is very broad with a maximum at oisyn ~ idt)~^ = j^uJc. Be- 
cause the synchrotron radiation is dominantly linearly po- 
larized, we consider E to be one of the linear polarization 
components. For the moment we neglect the influence of in- 
tervening plasma. The representation of E{t) by equation 
(1), where £j has a broad frequency spectrum, is specific to 
incoherent synchrotron radiation. It does not apply for ex- 
ample to thermal black-body radiation or inverse Compton 
radiation. 

We consider that the signal (1) is first passed through a 
bandpass filter which passes only frequencies in the interval 
/ ± _B/2, where B <C /. The resulting signal is then mixed 
with a precise local oscillator reference signal exp(27ri/f) to 
give 



£(f) = ^£:,A(f-i,) 



(2) 



where A(t) is the function which results from passing the 
delta function through the bandwidth B. For example, for 
a flat response within B, A(r) = sin(7r_Br)/(7rr) which has 
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a temporal width ~ . The expansion (2) is similar that 
of Appendix A of Cordes et al. (2004). 

The average power is proportional to an average of 
\E(t)\^ . We can assume {Sj£l) = {\£\'^)5jk because the dif- 
ferent electrons are uncorrelated. The individual pulses at 
times tj are closely spaced with a number per unit time J\f 
so that we can integrate over tj to obtain 

{\E\')=J^B{\£,f) . (3) 

We have used dt[A{t)f = B which holds for A(t) = 
sin(7rBf)/(7rt). 

Similarly, the autocorrelation functions is found to be 



- 1 r 



dtA{t)A{t + 5t) = pE{St) , (4) 



{E(t)E\t + St)) 

where the asterisk indicates the complex conjugate. Note 
that Pe{0) — 1 and that its width or correlation time is 
~ B^^. For the case A{t) = sm.{'nBt) / {-nt) one finds that 
PE{t) = A{t)/B. 

Consider the cross-correlation between the fields E\ and 
E2 at two separated frequencies /i and f-z ■ The signal E2 it) 
is obtained in a way similar to that for Ei(t), but it is mixed 
with a second local oscillator signal exp(27ri/2t), where the 
frequency /2 is accurately locked to /i. For simplicity we 
assume both bandwidths are equal to B. If the two frequen- 
cies are close together, (/2 — fif ^ /1/2, the amplitudes £j 
do not change significantly between the two frequencies be- 
cause of the broadband nature of the synchrotron radiation. 
Thus 



! PE{5i) , 



(5) 



{E^{t)E*2{t + 5t)) 

with pe given in equation (4). 

For the cross-correlation obtained by time-averaging 
over time intervals T, the variation in the difference between 
the two local oscillator frequencies, d(/2 — /i), must satisfy 

d(/2-/i)<(27rr)-i , (6) 

in order for the correlation not to bo smeared out. Thus the 
local oscillators will need to have very high stability possible 
with hydrogen maser sources. 

2.1 Intervening Plasma: Dispersion 

Propagation of the radio waves from a distant high redshift 
source such as a quasar will modify the cross-correlation be- 
cause of the frequency dependent delay time delay through 
the ionized intergalactic plasma. loka (2003) and Inoue 
(2004) discuss this intergalactic delay taking into account 
the redshifiting of the frequency and the time dilation dur- 
ing propagation through a homogeneously distributed fully 
ionized plasma assuming a fiat cosmology and express the 
delay as 



where DM is the dispersion measure. 



DM = 



dz'c 



dt 



dz' 



riejz') 
1 + 2' 



(7) 



(8) 



where \dt/dz\ = [{1 + z)H{z)\-'^ and H{z) = [0„(l-Fz)^-h 
Qa]^^^- Assuming the plasma is without loss or gain over 



time, the electron density varies as ne{z) = neoi^ + z)^ with 
rteo = const. Numerically, 



r(/)^4.14x 10-^':^ s 
Jio 



(9) 



with DM in units of pc(cm)~^ and /lo is the frequency in 
units of 10^" Hz. For redshifts z small compared with unity 

DM = f dine where dl is the path element to the source. 

loka (2003) and Inoue (2004) have evaluated the red- 
shift integral (7) for a homogenous intergalactic medium 
and estimate that DM « 1200« pc(cm)~^ for z < &. Addi- 
tionally, there may contributions to DM due to propagation 
through the free electrons of our Galaxy and/or the source 
galaxy. For our galaxy estimates of the DM in directions 
perpendicular /parallel to the disc plane range from 30 to 
lOOOpc(cm)"^ (Taylor et al. 1993). 

In order for the pulse-like nature of the synchrotron 
radiation not be smeared out by the variation of r(/) across 
the receiver bandwidth B, we must have B|dr/d/| < B~^ . 
This is the same as the requirement 



B < 1.1 X 10 



.3/2 

7 /lO 



Hz : 



.3/2 

3.16 X 10^^ Hz 



(10) 



DMi/2 

where the last expression uses DM w 1200a. We assume that 
B satisfies this condition. 

Consider again the cross-correlation of the fields at two 
nearby frequencies, /i and f2 = fi+ Sf, again with equal 
bandwidths B. We have in place of equation (5), 



{Ef{t-rf)E}+sf{t + 5t-rf+sf)} 



PEiSt-r'fSf) . (11) 



{\Ef\'} 

Here, Tfj+sf = T{f,f + Sf) and t} = dT{f)/df. The max- 
imum of pE occurs at 6t = r'fSf. Thus measurements of 
the values {5tm, Sfm) of the maximum allows one to derive 
Tf = 5tm/5fm which is directly related to the dispersion 
measure. 

Prom equation (9) we have 

-15 DM 



-8.28 X 10""°^ « -10"" -|- s^ 

JlO /lO 



(12) 



In order for the shift of the cross-correlation to be detectable 
with dispersion included we need |tj|(5/ > B~^ . In view of 
inequality (10) this implies that we need 



6f > 1.1 X 10^ 



j-3/2 
JlO 

DMi/2 



3.16 X 10 



,3/2 

,5 . 710 
.1/2 



Hz 



(13) 



which is the reverse of inequality (10). 

The small bandwidths indicated by inequality (10) may 
lead to poor signal to noise ratios in measurement of the 
cross-correlation functions. This may be offset by use of a 
filter bank with a large number iV ^ 1 of uniformly spaced 
channels each (j = L.A'^) of bandwidth B. The different 
channels can be averaged to give the cross-correlation as as 
'^j{Ej{t — Tf)Ej^f;{t + 5t — Tf+sf)) where j corresponds to / 
and k to 5f =const. Techniques and software for extracting 
weak dispersed pulses from radio measurements have been 
highly developed for applications to pulsar searches (Ran- 
som, Eikenberry & Middleditch 2002; Ransom 2012). 

A test and calibration of the proposed method can be 
obtained by applying it to the strong synchrotron emission of 
the Crab Nebula because in this case the dispersion measure 
is known from measurements of the Crab pulsar. It is DM 
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« 56.77 pc(cm) ^ (Lundgren et al. 1995). Inequality (7) 
then gives B < 1.46/fQ^MHz. The time delay in the cross- 
correlation between two frequencies separated by say Sf — 
5B is St = -3.44^is. 

2.2 Interstellar/Intergalactic Scintillations 

The interstellar scintillations (ISS) of small angular diame- 
ter extragalactic sources has been discussed by Lovelace & 
Backer (1972) and Goodman (1997). Lovelace and Backer 
considered the random modulation in the frequency domain 
of a wave from a broadband point source caused by the 
propagation through electron density irregularities in the 
interstellar medium (ISM). There is in general a "frequency 
correlation scale," A/, determined by the properties of the 
ISM. Two frequencies components of a wave from a point 
source separated by more than A/ are uncorrelated due to 
the random multi-path propagation. That is, the phase shift 
due to the extra propagation path length nfQgZ/c differs 
by more than 27r for frequencies separated by more than 
A/, where 6s is the scattering angle (proportional to /^^) 
and Z is the effective distance to the scattering screen. The 
value of A/ depends on the radio frequency and the propa- 
gation distance through the ISM. At frequencies / ^ /» ~ 2 
GHz, the scintillations are strong and A/// ~ (///*)^i 
whereas at higher frequencies the scintillations are weak and 
A/// ~ 1 (Lovelace & Backer 1972). Lovelace and Backer 
also estimated the critical angular size ^ of extragalactic ra- 
dio sources to be of the order of 25 x lO"'^ arc sec. Sources of 
smaller size may exhibit random intensity variations whereas 
larger sources are steady. The calculated shift in the cross 
correlation function arises from the contribution of the in- 
dividual radiating particles in the source and is not affected 
by the source size being larger or smaller than ^. 

Scintillations due to propagation through the inhomo- 
geneous intergalactic medium are considered by Goodman 
(1997) and estimated to be small compared with the ISS. 

In order for the shift in the cross-correlation function 
(equation 8) to be measurable we must have 5f < A/. This 
can be easily satisfied at GHz frequencies. 



3 CONCLUSIONS 

The detection of highly dispersed radio pulses likely of ex- 
tragalactic origin by Lorimer et al. (2007) and Keane et al. 
(2012) opens the intriguing possibility that the sources of 
such pulses can be identified and thus their distances de- 
termined. This would allow direct measurements of the dis- 
tributed baryon density of the intergalactic medium along 
lines-of-sight to possibly high redshift sources. This is of 
evident cosmological interest. Unfortunately, the narrow, 
highly dispersed radio pulses appear to be rare and non- 
repeating. Hence the distances to the sources are remain 
unknown. 

More recently, a search by Bannister et al. (2012) for 
prompt radio pulses at 1.4GHz from gamma ray burst 
sources - some at cosmological distances - has revealed 
highly dispersed pulses. However, the association of the 
pulses with the gamma ray burst sources remains uncertain. 

This work proposes a method to measure the disper- 
sion measures of high-redshift incoherent synchrotron radio 



sources associated with quasars by analysis of the amplitude 
cross-correlation of the radio waves received at two nearby 
frequencies. This method would allow direct measurements 
of the distributed baryon density of the intergalactic medium 
along lines-of-sight to sources of known high redshifts. This 
method utilizes the broadband nature of the synchrotron ra- 
diation. The bandwidths needed by this method are much 
narrower than those used to detect the mentioned dispersed 
radio pulses. Two key requirements on this method are iden- 
tified: one is that the local oscillators at the two frequencies 
have very high stability (inequality 6) and the second is that 
the bandwidths be sufficiently narrow (inequality 10). 
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